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We study reaction-diffusion processes on graphs through an extension of the standard reaction- 
diffusion equation starting from first principles. We focus on reaction spreading, i.e. on the time 
evolution of the reaction product, M(t). At variance with pure diffusive processes, characterized 
by the spectral dimension, d s , for reaction spreading the important quantity is found to be the 
connectivity dimension, d\. Numerical data, in agreement with analytical estimates based on the 
features of n independent random walkers on the graph, show that M(t) ~ t dl . In the case of Erdos- 
Renyi random graphs, the reaction-product is characterized by an exponential growth M(t) ~ e at 
with a proportional to ln(fc), where (k) is the average degree of the graph. 



A huge variety of different problems in chemistry, bi- 
ology and physics deal with reactive species in non triv- 
ial substrates pj. Seminal works on reaction and dif- 
fusion dynamics date back to the Fisher-Kolmogorov- 
Petrovskii-Piskunov (FKPP) model 



d t 9 = DA9 + f(0) , 



(1) 



where D is the molecular diffusivity, f{9) describes the 
reaction process and the scalar field represents the frac- 
tional concentration of the reaction products. Afterward, 
reaction-transport dynamics attracted a considerable in- 
terest for their relevance in a large number of chemical, 
biological and physical systems (jj . 

Complex networks are a recent branch of graph theory 
becoming very important for different disciplines ranging 
from physics to social science, from biology to computer 
science [3]. Although there exist an impressive amount 
of works on the study of both complex networks and 
reaction-transport processes, as far as we know, a gen- 
eral attempt to extend Eq. ([T]) on graphs and complex 
networks is still lacking. 

There are two main approaches to study reaction dy- 
namics on graphs. One concerns agent based models (La- 
grangian description) in which random walkers move on 
the graph and interact, with a given reaction rule, when 
they occupy the same site at the same instant [|| . A dif- 
ferent approach is based on a mesoscopic description of 
the reaction dynamics ((T|) in which diffusion is modified 
introducing a proper transport term taking into account 
the feature of the media in which the dynamic takes 
place Q. A particular approach in the mesoscopic de- 
scription of the dynamics (used, e.g., in the recent field 
of epidemic spreading Q), is to use a mean- field approx- 



imation in which a renormalized reaction term takes into 
account the network characteristics. 

The goal of this letter is to study the reaction spreading 
on graphs extending model ((T|). In the presence of more 
general transport processes, the diffusion term DAO in 
Eq. ([TJ, can be replaced by a suitable linear operator L. 
A general evolution equation for 9 is: 



d, 



l 



v, = L0+ -f(9) 

T 



(2) 



where we write explicitly the typical time scale, r, of the 
reaction process. An important class of processes of this 
type is the advection-reaction-diffusion (ARD), where 
L = -u-V+OA. Another interesting case is ruled by the 



effective diffusion operator L 



- 1 dr 



(k(r)i 



-1 d_ 
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suitable to study reaction dynamics on fractals [f|. 

Equation ([2]) is constituted by two terms: the trans- 
port term, L9, and the non-linear local reaction, f(9)/r. 
In the limit case without reaction, the link between the 
solution 0(x, f) and a suitable stochastic process is quite 
clear: for instance, if f(8) = and L = — u • V + DA, 
Eq. ((2|) is nothing but the Fokker-Plank equation asso- 
ciated to the Langevin equation dx/di = u + \j2Dr). In 
general, even in the presence of reactive terms and for 
general L, it is possible to write 0(x, t) in terms of tra- 
jectories using the Freidlin formula 



x, t) = ( 0(x, 0) exp 



1 f* f(9(x(s- 1 t),s)) i 
0(x(s;t),s) 



(3) 



where the average is performed over all the trajectories 
x(s;t) starting in x(0) and ending in x(t;t) — x. The 
possibility to write the generalization of ([3]) for a generic 
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diffusive process has been discussed in [9] . Following this 
approach we can determine the dynamical equation of 
reaction diffusion on graphs. 

As a diffusion process, we considered diffusion on an 
undirected, unweighted and connected graph G — (V, E), 
where V is the set of vertices of the graph (we consider 
a finite number, N, of vertices) and E is the set of edges 
connecting the vertices. The graph can be represented 
by its adjacency matrix A^ given by 10]: 



Aij — 



1 if (i,j)€E 
if {i,j)#E 



(4) 



The discrete Laplacian of the graph Ay [lOj, 111] is de- 
fined by: — Ay — fcj#y where fcj = J^-Ay, the 
number of nearest neighbors of i, is the degree of vertex 
i. Once the rate of the jump process, w, is introduced, 
the diffusion term can be written as 



dt 



(5) 



where 9i is the concentration at vertex i. Our goal is to 
add to this equation a reaction term: 



d0i A 



(6) 



The discrete-time version of the diffusion equation ([5]) is 
nothing but the random walk process described by the 
master equation 



(J) 



where jumps occur at time At, 2At, . . . , and the proba- 
bility for a walker being at vertex i to jump to the vertex 
j are given in terms of the adjacency matrix 



= w Aij At 



P 



% — y 1 

(At) 



if i^j 
1 — kiti) At if i = j 



(8) 



The discrete-time version of the reaction term can be 
defined as a non zero function only at discrete time when 
(5-form impulses occur: 



f(6,t)= g(0)S(t~nAt)At 



(9) 



where g(6) is a suitable function. Such a choice for f(6) 
allows us for a rigorous treatment of the discretization 
of the reaction term. With the above assumption, for 
Eq. ^ we have 



two steps: diffusion and then reaction The shape 
of the reaction map C?At(0) depends on the underlying 
chemical model. For auto-catalytic pulled reactions (the 
FKPP class, where, e.g., g(6) = 6(1 — 0)/t), character- 
ized by an unstable fixed point in 6 = and a stable one 
in 9 = 1 (the scalar field 9 represents the fractional con- 
centration of the reaction products; 9 = 1 indicates the 
inert material, 9 = the fresh one, and < 9 < 1 means 
that fresh materials coexist with products) one can use 
G A t{9) =0 + 6(1- 6)At/r. In the following we shall 
consider this type of reaction. 

The most important topological features of a graph 
can be related to the spectral dimension, d s , and the con- 
nectivity dimension, di , (also called chemical dimension) . 
The former is related to diffusion processes on graphs and 
can be defined in terms of the return probability Pa at 
site i for a random walker by d s = lim^oo — 2 ln ^" t ^ , 
or equivalently in terms of the density of eigenvalues of 
the Laplacian operator 11( . The connectivity dimension 
measures the average number of vertices connected to 
a vertex in at most I link, as #(Z) ~ l dl . For graphs 
embedded in an Euclidean space also the fractal dimen- 
sion df [l2| should be considered, describing the scaling 
of the number of vertices in a sphere of radius r in the 
Euclidean space, as #(r) ~ r df . The connectivity and 
fractal dimension can be different and they are related 
via the mapping between the two distances r and I [l3|. 

As a typical example of undirected, unweighted and 
connected graph, we show in Fig. [1] the reaction spread- 
ing in the T-graph [14| . The field 6 is initialized to zero 




FIG. 1. The percentage of quantity of product times r, M(t)r 
vs t. Numerical results for Eq. ([6]) with w = 0.5 are compared 
to prediction t '. For this graph d; = In 3/ In 2 ~ 1.585, di = 
2 In 3/ In 5 ~ 1.365. Inset: Spreading on a T-fractal. 



,(i + At)=G A t DOi(i) 



(10) 



where Ga*(#) = 6 + g(8)At is the assigned reaction map. 
It is worth noting that Eq. (|10[) can be seen as a numer- 
ical method when the integration of ^ is performed in 



in each vertex except the central one in which 8i (0) = 1 . 
Using Eq. ^ we study the time evolution of the system. 
An interesting observable to characterize the spreading 
of the reaction is the percentage of the total quantity of 
the reaction product, i.e., M(t) = jjJ2iev ^»'(*) wnere 
N is the total number of vertices. As clearly shown in 



3 



Fig. [H M(t) grows as a power law that can be interpreted 
as follows. Starting from a single vertex with 0j(O) = 1, 
after t step the number of vertices reached by the field 
is #(t) ~ t dl . Therefore, in the limit of very fast reac- 
tion, when each vertex reached by the field is immediately 
burnt (i.e, 9j — > 1), we can expect: 



M(t) ~ t d 



(11) 



Fig. [T] confirms that the connectivity dimension is the rel- 
evant quantity for the reaction spreading on graph. This 
behavior can be also understood thinking of the asymp- 
totic behavior of the reaction process as determined by 
the spreading of the front in the topological metric of the 
graph. In this case the characteristic time of reaction, r, 
appears only in the prefactor of the exponential. 

Moreover, a theoretical argument further confirms the 
importance of di. The analysis is based on an analogy 
between reaction spreading and short time regime of the 
number of distinct sites visited by n independent random 
walkers after t steps on a graph, S n (t) [T^]. This quan- 
tity can be computed as S n (t) = Y^jLo 1 — Coj(t) n , where 



Coj(t) is the probability that a walker starting from site 
has not visited site j at time t, the sum is over all the N 
sites of the graph dropping the dependence on the start- 
ing site 0. When the number of walkers is large (n — > oo), 
Coj(t) n tends to zero if site j has a non zero probability of 
being reached in t steps. In this limit, S n (t) represents all 
the sites which have nonzero probability of being visited 
by step t and, as t is equal to the connectivity distance, 
S n (t) ~ t dl . This is precisely the regime observed in the 
reaction spreading (see Eq. (fTTj) and Fig. [TJ. An esti- 
mate of the validity of the short time regime is given in 
terms of the smallest non zero occupation probability on 
the graph at time t, P m — (fc) _t , being (k) the average 
degree of the graph, i.e., the mean number of link for 
each vertex. As the short time regime is supposed to 
hold as long as nP m 1, one obtains that the reaction 
spreading regime is observed up to times t ~ Inn. On 
the other hand, the asymptotic regime is dominated by 
the number of distinct visited sites by a walker, that is 
S n (t) ~ t ds / 2 in graphs with compact exploration d s < 2, 
or simply by t on graphs with d s > 2 [15j | . In the case 
of very fast reaction regime, the front can be considered 
as equivalent to an infinite number of walkers, hence the 
asymptotic regime is never reached, leaving the dynamics 
to be governed by the sole di . 

As for the spectral dimension, it is the relevant quan- 
tity when dealing with random- walk dynamics Q and in 
some reaction diffusion processes. For instance in [l6j 
Eq. ^ has been studied for coarsening processes where 
f(0), at variance with our case, has a bistable structure. 
Moreover a further argument confirms the minor role of 
the spectral dimension in the case of reaction-diffusion 
dynamics using FKPP reaction terms. When dealing 
with standard diffusion ((x 2 (t)) ~ t) it is possible to 
show 9] that the spreading dynamics is the same dis- 
played by the standard reaction/diffusion problem (p}, 
i.e., M(t) ~ t d (where d is the dimension of the space). 



On the other hand, the presence of anomalous diffusion 
((x 2 (t)) ~ t 2v with v 7^ 1/2) does not implies that the 
spreading is anomalous: case exists Q in which diffusion 
is anomalous but reaction spreading is standard. 

The same behavior displayed in Figure [1] has been ob- 
served in several other self-similar graphs (e.g., Vicsek 
and Sierpinski carpet, not shown here), confirming the 
leading role of di. In the case of percolation clusters, the 
importance of the connectivity dimension, and the differ- 
ence between connectivity dimension and fractal dimen- 
sion was previously shown (l3j . 
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FIG. 2. Reaction product M(t) vs t for Erdos-Renyi graphs. 
Three results for different average degree of connectivity are 
shown in the fast reaction regime (r = 0.1). The reaction 
spreading follows an exponential behavior M(t) ~ t a , where 
a depends on (k) as shown in the inset. 



Now we focus on the behavior of Eq. $Q for Erdos- 
Renyi (ER) random graphs [l(| for which di — oo. In the 
ER graphs two vertices are connected with probability p. 
We choose p > so that the graph contains a global 

connected component. The average degree of the graph is 
(k) = p(N — 1) . On ER graphs the number of points in a 



sphere of radius / grows exponentially, 



hence 



we expect a similar behavior for the spreading process: 



M(t) 



(12) 



as shown in Fig [2j If (k) is large and the reaction is slow 
enough we have a two steps mechanism: first there is 
a rapid diffusion on the whole graph, then the reaction 
induces an increase of 9i. This leads to a simple mean 
field reaction dynamics, d t p(t) = p(t)(l — p(t))/r, where 
p is the average value of di on the graph. In this case 
a = 1/t as clearly observed in numerical simulations 
(not shown here). 

In the much more interesting case of fast reaction, at 
each time step the number of sites invaded is proportional 
to the average degree of the graph, so that after t steps 
we have: 



M(i) ~ (Ci(fc))* = e C3ln < fe > t ! 



(13) 
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leading toa~ ln(fc), see inset of Fig [2] 
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FIG. 3. The scaling exponent a normalized with ln(fc), as a 
function of the inverse of reaction-time -. The straight line 
indicate t 13 with 8 = -0.8. 



Furthermore, at variance with the case of graphs with 
finite di, at least in the case of fast reaction and FKPP 
reaction term, r plays an important role since C'2 is a 
function of r. Fig. [3] shows the dependence of the expo- 
nential behavior of reaction spreading rescaled with ln(fe) 
as a function of r. We can fit the dependence of the curve 



on r with a(k, r) ~ Ct 13 ln(fc), with /3 ~ -0.8. This scal- 
ing can be related to a mean field-like equation of the 
type: 



(14) 



We have considered a general model for reaction- 
diffusion dynamics on graphs, allowing for a general and 
detailed treatment of the diffusive and reaction terms. 
We study only large systems in which the asymptotic 
scaling for the reaction spreading is well defined. On the 
other hand, although the spreading dynamics on small 
systems is certainly a very interesting issue, it deserves 
careful attention and it is beyond the scope of the present 
work. In fact, even in the absence of reaction (i.e., pure 
diffusion) in small systems the boundaries can induce 
rather complicated behaviours [l7j . 

On undirected and finite dimensional graphs, we found 
that a major role in the reaction spreading is played by 
the connectivity dimension, which rules the asymptotic 
of the reaction product as a function of time. On random 
graphs with infinite connectivity dimension, the reaction 
spreading shows an exponential behavior, whose scaling 
depends on the average degree of the graph. In this case, 
we obtain two mean-field like equations, one in the slow 
reaction limit and one in the fast reaction limit. In par- 
ticular, in the fast reaction case, non-trivial dependence 
on both the average degree of the graph and the reac- 
tion characteristic time is shown. Our approach could be 
therefore suitable for a rigorous derivation of mean field 
like equations in more complex topologies. 
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